Breaking parameter degeneracy in interacting dark energy models from observations 
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We study the interacting dark energy model with time varying dark energy equation of state. 
We examine the stability in the perturbation formalism and the degeneracy among the coupling 
between dark sectors, the time-dependent dark energy equation of state and dark matter abundance 
in the cosmic microwave background radiation. Further we discuss the possible ways to break such 
degeneracy by doing global fitting using the latest observational data and we get a tight constraint 
on the interaction between dark sectors. 
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I. INTRODUCTION 

It is convincing that cold dark matter (DM) and dark 
energy (DE) are two dominant sources governing the evo- 
lution of our universe. Considering that these two com- 
ponents are independent is a very specific assumption, 
in the framework of field theory it is more natural to 
consider the inevitable interaction between them |l[ . An 
appropriate interaction in the dark sectors can alleviate 
the coincidence problem and understand why energy den- 
sities of dark sectors are of the same order of magnitude 
today 0-0. 

However, it was suspected that the coupling between 
dark sectors might lead to the curvature perturbation in- 
stability [7] . This problem was later clarified in Q where 
it was shown that the stability of curvature perturbation 
holds by appropriately choosing the forms of the inter- 
action and the values of the constant equation of state 
(EOS) of the DE. Another possibility to cure the insta- 
bility was suggested in 

The observational signature of the interaction between 
dark sectors has been widely discussed ^ was 

found that the coupling in dark sectors can affect sig- 
nificantly the expansion history of the universe and the 
growth history of cosmological structures. A number of 
studies have been devoted to grasp the signature of the 
dark sectors mutual interaction from the probes of the 
cosmic expansion history by using the WMAP, SNIa, 
BAO and SDSS data etc [lj]-[22j]. Interestingly it was 
disclosed that the late ISW effect has the unique ability 
to provide insight into the coupling between dark sectors 
21] . Furthermore, complementary probes of the coupling 
within dark sectors have been carried out in the study of 
the growth of cosmic structure [24[ -[111- ^ was found 
that a non-zero interaction between dark sectors leaves a 
clear change in the growth index [l5|, HU, [25| . In addition, 
it was argued that the dynamical equilibrium of collapsed 
structures such as clusters would acquire a modification 
due to the coupling between DE and DM [28|, [2!| , which 
could leave the signature in the virial masses of clusters 
[29| . The imprint of dark sectors interaction in the cluster 



number counts was disclosed in [271 ] . N-body simulations 
of structure formation in the context of interacting DE 
models were studied in 30]. 



Since both DE and DM are currently only detected via 
their gravitational effects and any change in the DE den- 
sity is conventionally attributed to its equation of state iu, 
there is an inevitable degeneracy while extracting the sig- 
nature of the interaction between dark sectors and other 
cosmological parameters. The degeneracies among the 
dark sector coupling, the the equation of state (EoS) of 
DE and the DM abundance were first discussed in [3l|. 
In the formalism of the perturbation theory, it was found 
that the degeneracy can be broken and tighter constraint 
on the interaction between dark sectors can be obtained 
from observations. 



In (3l|, the EoS of DE was taken to be constant. Re- 
cently, more accurate data analysis tells us that the time 
varying DE gives a better fit than a cosmological constant 
[32| . The time varying EoS influences a lot on the uni- 
verse evolution, perturbation stability and it must affect 
the constraint on the interaction between dark sectors. 
Thus it is of great interest to generalize the discussion 
in 31] to the interacting DE model with time-dependent 
DE EoS. The observational constraints on an interacting 
DE with time varying EoS were discussed in [19L |33| . In 
this work we will first discuss the stability of the pertur- 
bation and the degeneracy between the dark sector inter- 
action and the time varying EoS of DE. Furthermore we 
will explore possible ways to break the degeneracy among 
the dark sectors' coupling, the time-dependent DE EoS 
and the DM abundance. In our study we will use the 
popular Chevallier-Polarski-Lindcr (CPL) parametriza- 
tion [34| to describe the time varying DE EoS, which 
expresses the EoS in terms of the scale factor in the form 
w(a) — awo + (l — a)w e . In the early time, a <C 1, w ~ w e . 
At the present time when a ~ 1, w ~ wq. 
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II. PERTURBATION FORMALISM AND ITS 
STABILITY 



In the spatially flat Friedmann-Robertson- 
Walker(FRW) background, if there is an interac- 
tion between DE and DM, neither of them can evolve 
independently. The (non) conservation equations are 
described by 



p' c + 3W p c = aQ c , 
p' d + 3H(l + w)p d = aQ d , 



(1) 

(2) 



where the subscript c represents DM and d stands for DE. 
Q\ is the term leading to energy transfer. Considering 
that there is only energy transfer between DE and DM, 
we have Q c = —Qd = Q- The sign of Q determines 
the direction of the energy transfer. For positive Q, the 
energy flows from DE to DM. For negative Q, the energy 
flow is reversed. In our following study, we adopt the 
phenomenological interaction in the form Q = 3^Hp c , 
where £ is the coupling constant. This merely repeats the 
analysis of [10], where the energy transfer was assumed 
to be proportional to the energy density, but here we 
have the energy density times the Hubble parameter. In 
[7] , it was indicated that qualitatively similar conclusions 
apply to both cases. The presence of DM energy density 
in the interaction was shown problematic in the stability 
of perturbation Q • It was argued that only for constant 
EoS of DE in the phantom region, this interaction form 
can be viable and effective to alleviate the coincidence 
problem [sj. It is of interest to explore this interaction 
form when the DE EoS is time-dependent with the CPL 
parametrization, especially when the DE EoS is in the 
quintessence region. 

As discussed in [3l[ , the coupling vector is specified in 
co-moving frame as 



Qc. 



(A) 



(^,o,o,of. 



The perturbed form of the zero component SQ®^ can 
be uniquely determined from the background energy- 
momentum transfer and has the form 

SQ° (X) = -*Q (A )/a + SQ (x) /a. 

While the spatial component of the perturbed energy- 
momentum transfer SQ^ can be set to zero since there 
is no non-gravitational interaction in the DE and DM 
coupled system and only the inertial drag effect appears 
in the system due to the stationar y e nergy transfer be- 
tween DE and DM as discussed in [311 ] . 

The perturbation equations for DM and DE have the 



forms |2J, |31j 

D' c = -kU c + ZH^a -3£$' 



U' 



-HU C + - 3H£U C 



D' d = -3H(C 2 - w)D d - 9H 2 (C 2 - Cl)^- 

+ {3w' - 9H(w - C 2 ){£r + 1 + w)}$ 
+ 3^r$' - 3H^ + 3H£r(D d - D c ) 
U d 



kU d 



(1 + w)k 

U' d = -H(l - 3w)U d - 3C 2 (£r + 1 + w)fc$ 



+ 3U{C 2 e -C 2 a )U d 



+ 3U{Cl-Cl)tr-^L 
1 + w 

+ kC 2 e D d + (1 + w)k^ + 3H£rlI d . 

r = pd Pd is the ratio of DM to DE. In the above, C 2 = 
w < 0. However, it is not clear what expression we should 
have for C 2 . In Q it has been argued in favor of C 2 = 1 . 
This is correct for the scalar field, but it is not obvious 
for other cases. The most dangerous possibility, as far as 
instabilities are concerned, is C 2 — 1 ^ C 2 = w < since 
the first term in the second line of the equation of U d can 
lead to the blow up when w is close to —1. Assuming 
that C 2 — l,C 2 — w, the last two equations above for 
DE can be rewritten as 

D' d = 3U(-1 + w + $r)D d - 9H 2 (1 - w)(l + -£—)^r 

1 + w k 

- kU d + {3w' + 9U{1 - w){£_r + 1 + ?«)}$ + 3£r$' 

- 3H^ - 3H£rD m , 



U', 



-2H{1 



3£r 



1 



)U d + kD d - 3(£r + 1 + w )k$ 



(1 + w)k^f. 



Numerically, we find that the first two terms on the RHS 
of the perturbation equations for DE contribute more 
than other terms to the divergence. In order to explain 
the reason for the blow-up, we keep the leading terms 
and use ~ -w in the early universe @. In the early 
universe w ~ w e and the above equations can be reduced 
in the form 



D' d = -3UD d - 9H' 



1 - W e U d 



1 



k 



U' d = 2H \ 2We U d + kD d . 

1 + We 



The second order differential equation for D d can be writ- 
ten as 



n" 



In the radiation dominated period, we have H 
1/t,W 1/t 2 , thus 

l + 3w e D' d 6 
D d - ■ ^L) d , 

1 + W e T T z 
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which gives the solution 

D d = C lT ^ +C 2 T r ~, 

where 



r± = 



1 + Aw e + ±a/-5 - 4u> e + lOw? 



1 



When r± is a real positive value, the perturbation will 
blow up. In 8] for the constant DE EoS, w was required 
to be smaller than — 1 to accommodate the stability. For 
the DE EoS with CPL parametrization, in the early time 
w ~ w e < — 1 can also allow negative r± to provide sta- 
bility. Moreover as shown in Figfl] when w ~ w e > — 1/4 
in the early time, the stability can be safely protected as 
well. 
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The Behavior 
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meet at the point where the square root becomes zero. 
The real part of r± are the same on the right side of 
this point. 



When —1 < w e < —1/4, the above approximate an- 
alytical analysis cannot help, we need to count on the 
numerical calculation to see the stability of perturba- 
tion in the early universe. Results are shown in FigJ5] 
When w e falls in this range, the perturbation grows 
slower for weaker coupling. When the existing dark sec- 
tor coupling is weak enough, it has the possibility to keep 
the perturbation stable. For smaller w e in the range 
— 1 < w e < —1/4, weaker interaction is required to keep 
the stability in perturbation. 

Thus unlike the constant DE with w > — 1, the time 
varying DE can allow stable perturbation when it is in- 
teracting with DM. However in some ranges of w in the 
quintessence region, in order to keep the stability in the 
perturbation we have the upper limit on the strength of 
the interaction. When the DE EoS is in the phantom 
region, the stability of the perturbation is always pro- 
tected in the presence of the interaction between dark 
sectors, which is the same as the case we observed for 
the constant EoS of DE Q . 
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FIG. 2: The evolution of the curvature perturbation 
and DE density perturbation for various £ when 
-1 < w e < -1/4. 



III. DEGENERACIES OF INTERACTION, DE 
EOS AND DM ABUNDANCE IN CMB 
SPECTRUM 

With the perturbation formalism at hand, we can 
study the influence of the interaction between dark sec- 
tors and other cosmological parameters on the CMB 
power spectrum. We will concentrate our attention on 
the DE EoS in the quintessence region w > — 1, since 
it has different property as discussed in the above sec- 
tion compared with the case with constant EoS. In FigJ3] 
we illustrate the theoretical computation results of the 
CMB power spectrum for changing wo,w e in the CPL 
parametrization of DE EoS and coupling constant £ but 
with DM abundance fixed. 

We see in the CMB TT angular power spectrum that 
the change of the w e only modifies the acoustic peaks, 
while it does not influence the low-1 part of the spectrum. 
The influence of wo is on the contrary, it changes little 
of the acoustic peaks, but influences more on the small- 
1 spectrum. When wq decreases, low-1 spectrum will be 
suppressed. Different effects given by w e and wq on CMB 
spectrum are interesting. Since we have more CMB data 



4 



7000 
6000 
5000 
4000 
3000 



= -1-0. fl m 
; = 0.000, w e = 
; = 0.000, w a = 
; = 0.000, w e = 
: = 0.005, W. = 
; = 0.005, w s = 
; = 0.005, w e = 
; = 0.010, w e = 
: = 0.01 0, w = 
, = 0.010, Wf 



= 0.14 

0.0 — 

0.1 — 

0.2 --- 

0.0 - 

0.1 - 

0.2 ... 

0.0 - 

0.1 - 

0.2 --- 




(a) w e 




(b) w ~ § 

FIG. 3: The dependences of CMB angular power 
spectrum on the time varying DE EoS and the coupling 
between dark sectors. 



in large-1 than small-1 spectrum, it is easier to constrain 
w e than wo. This can be used to explain the fitting results 
in the following. 

The above discussion is for the fixed coupling constant. 
If we allow the change of £, we see that at the acoustic 
peaks, £ plays more important role in the change of the 
spectrum than the influence of w e . In the low-1 region, 
we also observe that £ effect is more important than wo 
to suppress the CMB spectrum. 

The above discussion is valid for fixed DM abun- 
dance. Now we investigate the dependence of CMB 
angular power spectrum on the abundance of matter, 
u! rn = fl m h 2 . Since the abundance of baryon is fixed, 
it is identical to investigate the effect of the abundance 
of CDM. Although the abundance of the DM does not 
affect much on the low-1 CMB power spectrum, it quite 
influences the amplitude of the first and second acoustic 
peaks in CMB TT angular power spectrum (see Figj4j) . 
Decreasing uj m will enhance the acoustic peaks. This ef- 
fect is degenerated with the influence given by the dark 
sectors' interaction and w e as we observed in Fig|31 A 
possible way to break this degeneracy is to consider the 
influence of the interaction on the low-1 CMB spectrum. 



Moreover, we can include further observations to get a 
complementary constraint on the matter abundance and 
this in turn can help to constrain the coupling between 
dark sectors. 




FIG. 4: The dependence of CMB angular power 
spectrum on the DM abundance. 



In order to extract the signature of the interaction 
and put constraints on other cosmological parameters, 
we need to use the latest CMB data together with other 
observational data. We report the results of fitting in the 
next section. 



IV. FITTING RESULTS 

In this section we confront our models with observa- 
tional data by implementing joint likelihood analysis. We 
take the parameter space as 

P = {h,uj b ,uj m ,T,ln[l0 lo A s ],n s ^,wo,w e ) 

where h is the hubble constant, ojf, = fl(,h 2 ,uj rn = Q m h 2 , 
A s is the amplitude of the primordial curvature pertur- 
bation, n s is the scalar spectral index, £ is the coupling 
constant when we choose the interaction in proportion 
to the energy density of DM. To avoid the negative en- 
ergy density of DE in the early time of the universe [23l|. 
we limit the coupling constant £ to be positive which 
indicates the energy flow from DE to DM. wq, w e are pa- 
rameters in the CPL parametrization of DE EoS. Since 
we concentrate on the DE EoS in the quintessence region, 
we set the limit —1 < Wq < 0, — 1 < w e < in our data 
fitting. We choose the flat universe with ilk — and our 
work is based on CMBEASY code^. The fitting results 
are listed in Table JT] 

When we only use the CMB anisotropy data from 
the seven-year Wilkinson Microwave Anisotropy Probe 
(WMAP), in Fig 15] we show the Id marginalized likeli- 
hoods for all of the primary MCMC parameters of our 
model. It is clear that, just from CMB data, the con- 
straint on wq is worse than that on w e . This is consistent 
with the analysis we did in the theoretical study. The w e 
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affects more on the acoustic peaks. More observational 
CMB data there can help to constrain w e well. For our 
model with the interaction between dark sectors in pro- 
portion to the energy density of DM, we find that CMB 
data alone can impose tight constraints on couplings and 
£l m h 2 . This can be understood by using our theoretical 
analysis that the degeneracy between the coupling and 
the DM abundance can be broken by looking at the low-1 
CMB spectrum. 

In order to get tighter constraint on Q, m h 2 , we use 
the BAO distance measurements [36[ which are obtained 
from analyzing clusters of galaxies and test a different 
region in the sky as compared to CMB. BAO measure- 
ments provide a robust constraint on the distance ratio 



r s (z d )/D v (z) 



(3) 



where DJz) = [(1 + z) 2 D 2 A z/ H(z)] 1 / 5 is, the effective dis- 
tance [37J, Da is the angular diameter distance, and H(z) 
is the Hubble parameter. r s {zd) is the comoving sound 
horizon at the baryon drag epoch where the baryons de- 
coupled from photons. We numerically find za using the 
condition f/R = 1,R 

Xbao 1S calculated as [3^ 



If-aa defined in |38|. The 



Xbao 



-* ->0bs rr , 

(d-d ) T C _1 (d-d 



obs „ 



(4) 



where d = (d z=0 . 2 , d z=0 . 35 ) T , d = (0.1905, 0.1097) 5 
and the inverse of covariance matrix read 36] 



30124 



-17227 



-17227 86977 



Furthermore, we add the BAO A parameter [SS 



(5) 



.4 



0.35 



0.35 



dz 



-i 2/3 



E(z) 

0.469(n s /0.98)" a35 ± 0.017 



(6) 



where E(z) — H Sr and n s are the scalar spectral index. 
In order to improve the constraints on the DE EoS w, 
we use the compilation of 397 Constitution samples from 
supernovae survey [40(. We compute 



XsN 



E 



(7) 



and marginalize the nuisance parameter. 
We implement the joint likelihood analysis, 



DE EoS are also improved. Compared with the WMAP 
data alone, we see that the joint analysis by including 
other observational data provides tighter constraints on 
the cosmological parameters. 

In Fig [6] we plot the 2d marginalized likelihood for the 
interacting models from our MCMC run with CMB and 
other observational data. The shaded regions are the re- 
sult from WMAP alone at 68% and 95% C.L. If we just 
look at the CMB result, we see that when w e < —1/4, 
the allowed £ becomes smaller when w e decreases. When 
w e > —1/4, the allowed £ is bigger. This fitting result 
supports our theoretical observation obtained in the sta- 
bility analysis. The dashed and solid lines are the 68% 
and 95% C. L. regions in combination of WMAP+SNIa 
and WMAP+BAO+SN data sets, respectively. Figf5] 
summarizes our findings for the possible relation between 
w e and the coupling constant from data fitting. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper we generalized our previous study [3l[ on 
the interacting DE model to the DE with time varying 
EoS. We reviewed the perturbation theory and examined 
the stability in different model parameters. Based upon 
the perturbation formalism we have studied the signa- 
ture of the interaction between dark sectors from CMB 
angular power spectrum. Theoretically we found that 
there are possible ways to break the degeneracies among 
the interaction, parameters in CPL paramctrization of 
DE EoS and DM abundance. This can help to get tight 
constraint on the interaction between DE and DM. 

We have performed the global fitting by using the CMB 
power spectrum data from WMAP7Y results together 
with the latest SNIa, and BAO data to constrain the in- 
teraction between DE and DM and other cosmological 
parameters. As anticipated from our theoretical anal- 
ysis, the global fitting can really break degeneracy of 
cosmological parameters and give tighter constraints on 
the interaction between dark sectors, DE EoS and DM 
abundance. The relation between w e and the coupling 
obtained from the fitting supports the result got in the 
stability analysis. 

With the successful experience to deal with the inter- 
acting DE model with time varying DE EoS, our next 
step is to extend our study to the field theory based 
model to describe the interaction between dark sectors. 
The preliminary attempt was carried out in l]. More 
efforts are required on this direction. 



X 



2 ,2,2 

Xwmap + Xsn ' Xbao- 



(8) 



The cosmological parameters are well constrained. When 
the coupling between dark sectors is proportional to the 
energy density of DM, its constraint is much tighter than 
that from CMB data alone. The likelihoods of the fitting 
results for the matter abundance, parameters WQ,w e for 
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FIG. 5: 1-D marginalized likelihoods for cosmological parameters. WMAP alone (red), WMAP+BAO(green), 
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